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1. Introduction 

Semi-inclusive deeply inelastic scattering (SIDIS) experiments are sensitive to many correla- 
tions between the direction of parton and hadron spins and intrinsic transverse momenta. These 
have been parameterized in a systematic manner using transverse momentum dependent parton 
distribution functions (TMDPDFs), see [I]. TMDPDFs describe the distribution of partons car- 
rying a longitudinal momentum fraction x and an intrinsic transverse momentum kj^ in a hadron 
as illustrated in Fig. 1. Here we give an update of our effort [2] to develop techniques suitable 
for the calculation of moments of TMDPDFs on the lattice. Note that TMDPDFs are not to be 
confused with generalized parton distribution functions (GPDs), which provide probability distri- 
butions with respect to the impact parameter rather than k\^. For an overview of recent hadron 
structure studies in lattice QCD, we refer to [3]. 

Fig. 2 illustrates the factorization of SIDIS into perturbative and non-perturbative parts. The 
lower blob represents the non-perturbative contribution of the nucleon and is described by 

^i^]{x,kr,P,s)^ydk-j^e-''' {p,s\ q{i)r^.^^f^o)q{o) \p,s), (i.i) 

where \P,S) represents a nucleon state of momentum P and spin S, F is a Dirac matrix and k is the 
quark momentum, with k^ = xP^ . The Wilson line '^^(^^o) connecting the quark operators ensures 
gauge invaiiance. In SIDIS, '^^(^^o) = ^[f/+~n_]'^[£+oo(7_,oon_]^[Don_,o] ^ concatenation of three 
straight Wilson lines running to light cone infinity and back [4, 5]. 




Figure 2: Factorized tree level diagram of SIDIS. 



Some examples of TMDPDFs are f\, gn and gir, defined via [1] 

^^^'^{x,kr,P,S) = Mx,kl) (1.2) 
^^y'^H^,kr,P,S) = ^5+gi.(x,^i) + '^g,r{x,kl). (1.3) 

2. TMDPDFs from the Euclidean lattice 

In this exploratory study, we simpUfy the calculation and employ a single straight Wilson line 
"^^(^,0) = '^[ifi] running from to i. In this case, the matrix element appearing in eq. (1.1) can be 
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evaluated directly on the Euclidean lattice, as long as we set f = 0. To facilitate the translation of 
the correlators evaluated on the lattice into TMDPDFs, we first parameterize the matrix element in 
terms of Lorentz-invariant amplitudes A,(£^, ^ • P), for example 



18 «M 
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(2.1) 
(2.2) 



(P, S\ q{i) 7^ ^ ^(0) \P, S)=4A2P^,+ AiniN^ A3 , 
{P,S\ q{t)y^r"i/ q{Q) \P,S) = -AniN A(,S^-AimN A-, P^l-S + AniN'^ A 

The amplitudes A, , l-P) are extracted on the lattice and then Fourier transformed into TMDPDFs. 
For example, from eqns. (1.2), (1.1) and (2.1) we get 

7o(\/^|^±|)2A2(£2,^P), 
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where Jq is a Bessel-function. The restriction to ^ = on the lattice translates into the constraints 



<0, \£-P\<\i\\P\ 



(2.4) 



which preclude us from evaluating the full x- and A:j^-dependence directly as in eq. (2.3), but are 
harmless if we are only interested in the first Mellin moment, i.e., if we integrate over x. For 
example, we obtain 



dxf\{x,k\) 
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(2.5) 
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3. Simulation technique and parameters 
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Figure 3 



The matrix element (P,S| q{t)T ^\^i^q{Q) \P,S) is evaluated us- 
ing ratios of three- and two-point functions as described in Ref. [2]. 
The non-local operator = q{i)r^[£Q] q{0) inserted in the three- 
point function contains the Wilson-line "^[cfi] , which is implemented 
as a product of link variables. For oblique angles, we approximate a 
straight line by a step-Uke path as illustrated in Fig. 3. 

For our studies we work with MILC gauge configurations [6] based on an AsqTad improved 
staggered quark action with 2-1-1 flavors on a 20^ x 64 lattice with a lattice spacing a « 0.12fm and 
a strange quark mass anis = 0.050. We have used three different light quark masses, awi„^ = 0.020 
(m^ w 500 Me V, 239 configurations), am„ j = 0.030 (mj^ ~ 600 MeV, 281 configurations) and the 
three-flavor degenerate case airiu^d = 0.050 (niji 760 MeV, 213 configurations). 

The gauge configurations have been HYP smeared and bisected in the temporal direction to 
double statistics. We are using domain wall propagators and sequential propagators previously 
calculated by the LHPC collaboration on these configurations, with the valence quark mass tuned to 
match the staggered sea (see, e.g., [7]). The sequential propagators feature a source-sink separation 
of 

4mk ~ ^source — 10, and are available for two lattice nucleon momenta P — (0,0,0) and P — 
(—1,0,0), the latter corresponding to 500 MeV in physical units. We neglect contributions from 
disconnected diagrams. We have developed our software using the Chi^oma library [8]. 
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4. Renormalization of the Wilson Lines 

The Wilson line ^ycx)^ in our non-local operator gives rise to a linear divergence, which has 
to be removed by a renormalization constant 5m proportional to the cutoff, given by a^^ on the 
lattice. Refs. [9, 10] show within continuum theory that the renormalized operator is of the form 

=Z-'exp(-5mL) (4.1) 

Here subsumes renormalization factors associated with divergences at the end points, and L is 
the total length of the smooth Wilson line. 

In lattice QCD, the linear divergence has been a long standing issue in the context of heavy 
quark propagators [1 1]. We have calculated aSm for link paths on the axes in leading order pertur- 
bation theory, adapting the procedure in Refs. [12, 13] to our action according to Refs. [14, 15]. 

However, perturbation theory is not expected to give accurate results. Therefore we have also 
sought to determine aSm non-perturbatively with our "taxi driver method", which is based on the 
assumption that on the lattice, L in eq. (4. 1 ) is given by the total number of link variables, and thus 
allows us to deduced the renormalization constants from the comparison of straight and step-Uke 
Unk paths. There are two variants of this method, one based on data from Wilson lines in a Landau 
gauge fixed ensemble, and one based on Wilson loops. For the moment, we assume that the quark 
mass dependence is weak, and use the constants determined on the three-flavor degenerate lattice 
also for the lighter quark masses. We are currently testing the validity of the taxi driver approach 
and alternative methods on several different lattice spacings. 

5. Results with Preliminary Renormalization 

Here we present some results obtained by applying the techniques sketched above. Note that at 
the present stage we regard our renormalization procedure still as preliminary. In the following, we 
label our distributions "sW" to indicate that they are based on straight Wilson lines and are therefore 
not strictly identical to the TMDPDFs defined and used in the literature and for the description of, 
e.g., SIDIS. 

In Fig. 4a we display results for 2A2{l'^,i-P = 0). They are obtained from a three-point func- 
tion with the operator (f), where the Wilson line is renormalized using different approaches 
as described above. The overall normalization can be obtained requiring charge conservation, 
242"''(0,0) = 1. Applying eq. (2.5) to the fit results, we get /^'''^(^t^^) as plotted in Fig. 5a, which 
is interpreted as the unpolarized distribution of quarks in the unpolarized nucleon in Fig. 6a. In 
the unpolarized channels, the nucleon looks axially symmetric. For the large pion masses currently 
analyzed, the quark mass dependence of the width of this distribution appears to be rather weak, 
see Fig. 5c. 

The axial symmetry is distorted in the polarized case. As an example, we can consider the 
distribution of quarks with positive helicity A=+l, con^esponding to an operator (ff- {€) with F = 
7+^(1 + 7^) in a nucleon transversely polarized with transverse spin S\_. Using eqns. (1.2), (1.3) 

Pik^s^) ^ C j.<i>[^^(>+^)/2i(-,^^;P,5) = \ + — g^KA)] ■ (5.1) 
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(a) (b) 




Figure 4: (a) Amplitude 2Re A2{i'^,i-P = 0) for up quarks. The continuous curves are fits of tlie 
form Ci exp{—£^/af) +C2exp(— ^^/a|). We liave used tliree different sets of renormal- 
ization constants: For the upper curve, we have determined a5m from 1-Ioop perturbation 
theory, for the two lower curves we have employed two variants of the taxi driver method, 
(b) The same for the amplitude -2ReAT{e^,i-P = 0). 



This density is axially asymmetric due to the contribution from the TMDPDF gyj, which contains 
information about the correlation of spins and momenta of the form s-P k-S X ^j^ .Sj^. According 
to eq. (2.6), we obtain g[j^^{k^) from the amplitude A-j, as shown in the case of up quarks in 
Fig. 4b and 5b. For down quarks, the results are of opposite sign and smaller by a factor of about 
1/6. In terms of the density p(^j^,5'j^), we find that quarks of a specific polarization can have a 
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Figure 5: (a) /{ (k^) for up quarks as obtained from the Fourier transform of the fits in Fig. 4a. 
mtr^iT^*(^i) ^'^^ quarks as obtained from the fits in Fig. 4b. 
(c) Linear chiral extrapolation of the root mean squared transverse momentum {k\)^^'^ 
from fl^^^^{k\) for up quarks minus down quarks. 
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(5.2) 



non- vanishing average transverse momentum 

/ cP-k±^ k±p{k±,S±) 
Jd^k^p{k^,S^) 

and that it is opposite in sign for up- and down quarks, see Figs. 6b and 6c. That such deformed 
quark densities are to be expected has been mentioned, e.g., in Ref. [16]. Similar deformations 
have also been observed in the framework of GPDs in [17, 18]. 



6. Testing factorization in x and k± 

In previous sections we have studied lat- 
tice data for £-P = 0. Let us now explore the 
£-P-dependence of A2(^^,^-P). The first in- 
tegral in eq. (2.3) shows that it is related to 
the Bjorken-;c-dependence of ff^{x,k±) via 
a Fourier transformation. Figure 7 gives an 
overview of the unrenormalized data avail- 
able for Re A2 for up minus down quarks. 
The sector with data points is constrained by 
eq. (2.4) and the largest available nucleon mo- 
mentum \P\ of about 500 MeV. Figure 8 dis- 
plays the ^-P-dependence of the real and imag- Figure 7 
inary parts of A2 and reveals consistency of 

our results with Ai{i^,£-P) = (£, —l-P), which follows from the transformation property of the 
matrix element under Hermitian conjugation. 

In phenomenological applications, it is often assumed that f\{x,k\^) factorizes into an x- 
and a -dependent part, see, e.g., [19]. In our case, the hypothesis that ff^ is of the form 
ff{x,ki_) = lf{x) fi^'^{ki_) translates mloA2{f,£-P) = A2(£-P) A2{f,0) using the Fourier 
transform eq. (2.3). To test this hypothesis, we introduce a scaled amplitude 




(6.1) 



ReA2(£2,0)' 

Note that renormalization factors cancel in this ratio. If factorization holds, A2{i^,i-P) should 
be £^ -independent. We plot this quantity in Fig. 9a, selecting the imaginary part as an example. 
No significant ^^-dependence is visible, i.e., we confirm the factorization hypothesis within the 
accessible kinematic range and within our statistics. Given that A2{i^,i-P) is approximately 
independent, we can plot it with respect to £-P as in Figs. 9b and 9c. It is interesting to see that the 
result is qualitatively similar to the Fourier transform / dxexp{i£-Px)fi (x) of a phenomenological 
parametrization of the PDF /i (x), such as the one provided by CTEQ5 [20]. 

7. Conclusions and Outlook 



Presently, we employ a straight Wilson line between the quark fields in the definition of TMD- 
PDFs. We have shown first preliminary results from lattice QCD for the TMDPDFs ff^ and gf^ as 
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Figure 6: Quark density plots. Here the renormalization constants have been chosen according to 
the taxi driver method based on Wilson lines. 

(a) = / dx^^'^^\x,ki^ for up quarks at mj^; = 500 MeV. We interpret this as 
the charge density of up quai^ks in the nucleon in the transverse momentum plane. 

(b) Density p{kj^,Sj^) of up quarks with positive helicity A = + 1 (i.e., with spin point- 
ing in z-direction) in a nucleon polarized in transverse x-direction = ( 1 , 0) , evaluated 
at a pion mass niji = 500 MeV. The distribution features an average transverse momen- 
tum shift (^v) = (67 ± 5stat i 3ienorm.) MeV, whcrc the uncertainty from renormalization 
has been estimated from the comparison of the three different sets of renormalization 
constants employed in Figs. 4 and 5. 

(c) Same as in (b) but for down quarks. The average transverse momentum shift (^v) = 
(—24 lb 5stat i 3renorm.) McV has the oppositc sign as for up quarks. 

(d) Same as in (b) for up quarks minus down quarks. The deformation appears amplified. 
Note that this is not a density and not necessarily positive. 



a function of transverse momentum. We find that densities of longitudinally polarized quarks in a 
transversely polarized proton are deformed. Moreover, we confirm that the factorization hypothesis 
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ff^{x,kj^) = fj*(;c) f\^^^'^{kj^) is valid within the statistics of our data set. 

Concerning our renormaUzation procedure, further investigations are in progress. Further- 
more, we would like to extend our work towards non-straight Wilson lines, similar to those appear- 
ing in the definition of TMDPDFs for experimental processes such as SIDIS. 
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